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Abstract 

The quasiclassical Green functions of the Dirac and Klein-Gordon equations in 
the external electric field are obtained with the first correction taken into account. 
The relevant potential is assumed to be localized, while its spherical symmetry is 
not required. Using these Green functions, the corresponding wave functions are 
found in the approximation similar to the Furry-Sommerfeld-Maue approximation. 
It is shown that the quasiclassical Green function does not coincide with the Green 
function obtained in the eikonal approximation and has a wider region of applicability. 
It is illustrated by the calculation of the small-angle scattering amplitude for a charged 
particle and the forward photon scattering amplitude. For charged particles, the first 
correction to the scattering amplitude in the non-spherically symmetric potential is 
found. This correction is proportional to the scattering angle. The real part of the 
amplitude of forward photon scattering in a screened Coulomb potential is obtained. 



1 Introduction 



As known, the use of the wave functions and Green functions of the Dirac equation in an 
external field is a convenient tool for the calculation of QED amplitudes in the field. At 
high energy and small scattering angles, the characteristic angular momenta of involved 
particles are large, and the quasiclassical approximation is applicable. In this case the use 
of quasiclassical Green functions greatly simplifies the calculations. 

The quasiclassical Green function G(r 2 , r^e) of the Dirac equation in a Coulomb poten- 
tial was first derived in |2|] starting from the exact Green function of the Dirac equation 
||. The more convenient representation of the quasiclassical Coulomb Green function for 
the case of almost collinear vectors ri and r 2 was obtained in M, [jj . In the same geometry 
the quasiclassical Green function in arbitrary spherically symmetric decreasing potential 
was found in || [/J . 

In section 2 we derive the quasiclassical Green functions of the Dirac and Klein-Gordon 
equations for the case of a localized potential not assumed to be spherically-symmetric. 
We use the term 'localized potential' for the decreasing potential having a maximum at 
some point. The Green functions are obtained with the first correction taken into account. 
With their help the quasiclassical wave functions of the Dirac and Klein-Gordon equations 
and corrections to them are found (section 3). These wave functions generalize the results 
obtained in the Furry-Sommerfeld-Maue approximation || || pTO| . 

In the calculation of amplitudes of the high-energy processes the eikonal approxima- 
tion is often used. The corresponding wave functions and Green functions differ, generally 
speaking, from the quasiclassical ones and have a narrower region of applicability. There- 
fore, the use of the eikonal approximation without a special consideration of its applicability 
may lead to incorrect results. As an example, we show this in section 4 for the problem of 
small-angle scattering of a charged particle in an external field. In this section we present 
also the consequent derivation of the expression for the scattering amplitude using the qua- 
siclassical wave function. In particular case of a Coulomb potential the use of the eikonal 
wave function instead of the quasiclassical one leads to the incorrect result. 

In section 5 the quasiclassical Green function obtained is used in the calculation of 
the amplitude of forward elastic scattering of a photon in the atomic electric field (forward 
Delbruck scattering). As shown in 0, at this calculation it is necessary to take into account 
the correction to the quasiclassical Green function. This correction should be taken into 
account in the region, where the eikonal approximation is valid. The contribution of higher 
orders of the perturbation theory with respect to the external field (Coulomb corrections) 
is determined by the region, where the eikonal Green function is inapplicable. The real part 
of the forward Delbruck scattering amplitude for a screened Coulomb potential is obtained. 
It becomes comparable with the amplitude of Compton scattering already at relatively 
small energies of a photon. It may be important for the description of photon propagation 
in matter. 
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2 Green function 



As shown in papers |], [|, in the calculation of amplitudes of various QED processes it 
is convenient to use the Green function of squared Dirac equation D(r 2 , ri|e), which is 
connected to conventional Green function G(r 2 ,ri|£) by the relation 

£(r 2 ,n|e), (1) 

where 7 M are the Dirac matrices, p = — zV is the momentum operator, and V(r) is the 
potential. In the quasiclassical approximation it is possible to present function D(r 2 ,r 1 \e) 

as 
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(3) 



and k 2 = e 2 —m 2 . Thus, the problem is reduced to the calculation of the quasiclassical Green 
function of the Klein-Gordon equation with potential V(r) (Schrodinger equation with 
Hamiltonian H). 



Let us pass in function D^(r 2 ,r 1 \ e) from variables 17 and r 2 to variables 17 and r 
r 2 — 17. In terms of these variables the function D satisfies to the equation 



2/#(i7 + r) 
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ri,ri e) 
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where 
and A 



XV — V 2 /2k, X = e/k, p = — zV r . In the ultrarelativistic case A = +1 for e > 
-1 for e < 0. 



We seek for a solution to this equation in the form 

n(nl , , N exp(mr) 



47rr 



F(r,n) . 



(5) 



Note, that the factor in front of F in (^|) is a Green function of the equation (Q) for = 0. 
The function F satisfies the equation 



2k 



rA(F/r) 



(6) 



with the boundary condition F(r = 0, 17) = 1. For further needs we define an effective 
impact parameter of a rectilinear trajectory I, connecting point 17 and r 2 , as 



p* 



mm 



(7) 



Let this minimum is achieved at some point x = r*. The necessary condition Kp* ^> 1 for 
the applicability of the quasiclassical approximation is assumed to be fulfilled. Introducing 
the notation 07 2 = |i7 2 — r*|, we have two overlapping regions: 



1) min(a 1? a 2 ) <C Kp\ 

2) min(ai, a 2 ) » p* . 



(8) 
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In the first region, the r.h.s. of the equation (||) is small. The solution to this equation 
with zero r.h.s. has the form 



F = exp 



-zr 



/ 4>{ri + xr)dx 
Jo 



(9) 



which corresponds to the eikonal approximation. For the calculation of the first correction 
to F we search for a solution of the equation (|) in the form F = F (l + g) and neglect g 
in r.h.s.. As a result, we obtain the following equation for g: 



2in-^-g = (f) 2 (ri + r) + ir J dxx 2 A 1 (f)(ri + xr) + r 2 J dx xVi_i_0(r! + xr) 



(10) 



where the index 1 at the derivative designates differentiation over ri, Vu_ is a component 
of a gradient, perpendicular to r, Ai = V 2 . Integrating over r we obtain for g 



2k 



r 2 J dx x(l — x)A 1 (f)(r 1 + xr) — ir J dx 4> 2 {r\ + xr) — 
. o o 

1 x 

— 2ir 3 J dx(l — x)Vu_0(ri + xr) J dy y'V i±(j)(ri + yr) 
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Finally, up to the first correction, the Green function in the first region reads 



(11) 
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%kt — iXr / V(ri + xr)dx x 
l 

Ar 2 y dxx(l — x)AiV(r 1 + xr) — 



(12) 



- 2ir 3 J dx{l - x)Vi_L^(ri + xr) J dyyV^V {ri + yr) 



In this expression we put A 2 = 1 assuming e ^> m. Substituting (^3|) in (||), we find the 
expression for function D in the first region 



■D(r 2 ,ri|e) 
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exp 
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x<l — — [ dxcx ViV(ri + xr) + ^— [ dxx(l — x)AiV(r a + xr) — 
2k J 2k J 

K o o 

• 3 1 x 

— — J dx(l — x)Vi_Ly(r! + xr) J dyy'Vi ± V(r 1 + yr) 



In the second region the r.h.s. of the equation 
nates, we can rewrite (^[) as 



is not small. Using spherical coordi- 



.d_ 

Or 



- 0(ri + r - 
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where L 2 is the angular momentum operator squared. In this equation r is a free variable. 
We are interested in a value of function F(r, i^) at r = r 2 — r x . It is convenient to direct the 
polar axis along r 2 — ri. The term in the l.h.s. of flTil) containing L 2 should be taken into 
account for r > ai. It can be checked by applying the operator L 2 to the eikonal function 
(§). When calculating the function -F(r 2 — ri, ri) it is enough to consider a narrow region 
of polar angles of vector r where 9 ~ p*/ai <C 1. The r.h.s. of (14) is small. We seek for a 
solution to this equation in the form 



1 1\ L 2 

r ai/ 2k 



(15) 



Substituting flTJI) in (|l4l), we obtain J 7 



<9r 



T=- — 
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where 
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.-iA, 



zL 2 fd_ 

Kr 2 \ Q r 



iL 2 ' 
4/tr 2 
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(16) 



(17) 



We are going to solve this equation up to the first correction in the parameter For 
this purpose we should keep only two terms in the expansion of in terms of commutators 
of the operator A and <fi: (j) ~ — i[A, <f>]. In the first approximation it is possible to 
neglect r.h.s. of fll6|) and to replace <ft by 0. Then the function T coincides with the eikonal 
function F (see (|9])). To find the first correction, we present T as T = F (l + g{). We 
obtain the following equation for gi 



n . d 
2%K—g x 
or 



n+r - -- 
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L 

Integrating it over r, we find 

i i 
J dx(p 2 (ri + xr) + ir J dxx (^1 x^j Aij_</>(ri + xr) + 

o fll 
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(19) 



Here = 'V 2 ± . With the help of the expansion in spherical functions it can be shown, 
that at (3 <C 1 for arbitrary function g(r) the following relation is true with the appropriate 
accuracy: 



exp[-i(3 2 L 2 } g(r) w / ^e iq2 g(r + 2(3rq 
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(20) 



where q is a two-dimensional vector, perpendicular to r. Using QlB]), ([19|) and (pOD, we 
obtain in the second region the following expression for the Green function with the 
first correction taken into account: 



D (0) (r 2l rik) 
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+ 2r 2 J dx (l - ^ar) (V 1± U(r x )) | dyyV 1± V(r v 



where r x = I - ! + xr + qxy2r(r — ai)/(reai), r = r 2 — ri. Let us remind, that the derivatives 
over ri in this formula must be calculated at fixed r. 

Note that the term proportional to q in r x is essential only in a narrow region \x—a\/r\ ~ 



p*/r <C 1. Using this fact, we can eliminate the quantity a\ from the formula (0). In order 
to do this, we present r x = H x + 5r x , where 



2r 



R-x = ri + xr + J 2x(l — x)r//tq, 5r x — \ — \x\jrja\ — 1 — J x(l — x) ) q . (22) 



K. 



expand V(r x ) in the exponent with respect to ^r^ up to the first term, and replace r x with 
R x in the rest of the expression. After the integration by parts over q we obtain: 
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(23) 







As mentioned above, the regions of applicability of two formulas (|T|) and (^) are over- 
lapping. At p* <C min(a 1 ,a 2 ) <C Kp 2 it is possible to expand V(R X ) in q up to the second 



term and integrate over q. After that the formula ( p3|) turns into ( [12|) including the terms 
of the order of 1/ftp*. Therefore, the result can be presented in the form which is correct 
everywhere for ftp* 1: 
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2 J dxV(K x ) - V( ri ) - V(r 2 ) 
o 

+ l -!ljdxjdy y x (i- x )y(i- y )-(i- x )y\ ( V 1± V (R x )) ( V 1± V(R y ) ) 



(24) 



x 1 - — 
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R x is defined in (p2|) . Indeed, in the first region we can expand V(Rx) in q up to the 
second-order terms and to integrate over q, which leads to the formula fljjD . In the second 
region one should bear in mind, that the main contribution to integrals comes from a 
narrow region of x and y close to a±/r with Sx ~ by ~ p*/r. Therefore 



yjx(l - x)y{\ - y) - (1 - x)y = x -y - (x -y) 2 - (j x{\ - x) - yjy(l - y)) 



x-y 



up to the second-order terms in p*/r. Besides, in this region the terms linear in V can be 
neglected in (| 



Thus, in the second region (]24|) turns into 
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Using (|2]), we obtain the final expression for the Green function -D(r 2 , rx|e) 



D(r 2 ,r 1 \e) 



4ir 2 r 



<iq exp 



iq 



+ 



iXr [ dxV(R x ) x (25) 
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1-L-Jdxa V 1 \/(R X ) - A 2 1 da^R,) - ^(n) - V(r 2 ) 
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2 1 x 

+ Jdxjdy y x (l-x)y(l-y) - (1 - ar)y] (V 1± U(R x .))(V 1± U(R y )) 

o o 

The advantage of the representations (]24]) and (p5|) is that they keep their form in any 
reference frame. 

The integration over the variable q in formulas (21) and (23) can be interpreted as the 
account for quantum fluctuations near the rectilinear trajectory between vectors ri and r 2 . 
The integral over q converges at q < 1. Using this fact, we conclude that the quantum 
fluctuations can be neglected, when 

J dx\J x{l — x)r/ k\Vh_V{yi + xr)\ <C J dx\V{ri + xr)\ 

This condition is actually equivalent to the first condition in (H) ensuring applicability of 
the eikonal approximation. 



In the formulas (|24|) and (|25|) the terms containing a potential in the pre-exponent 
factor, give the correction to the quasiclassical Green function. The expressions obtained 
are valid when these correction are small. Our results were obtained for a localized po- 
tential. Nevertheless, they are valid for any potential if the correction is small, e.g. for a 
superposition of localized potentials. 



2.1 Quasiclassical Green function in a central field 

For a spherically symmetric potential the quasiclassical Green function without corrections 
can be obtained also from the results of || 0, where it was calculated for the case of 
almost collinear vectors ri and r 2 . In these papers with the help of quasiclassical radial 
wave functions the following expression for the function D^°- , (r 2 , ri| e) was obtained at small 
angle 9 between vectors r 2 and — ri: 

D^(r 2 ,r 1 \e) = ^ dllJ (W) exp i \ 1 + + 2XS(1/k) + A(<fr(ri) + ^(^2)) • 



(26) 

Here 



$(r) = / V(()d( , 5(p) = ~Jv + <K, A = e/ K . 

r 

If the angle 6' = 71 — 9 between vectors ri and r 2 is small, then 

1 

47r|ri - r 2 | 



£> (0) (ri, r 2 | e) = - f exp{m|ri - r 2 | + iAsign(ri - r 2 )($(ri) - <&(r 2 ))} . (27) 
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For further transformations it is convenient to rewrite the expression 
using the identity 

dllJ (W)g(l 2 ) = i-| dqexp^qfl)^ 2 ), 



in the other form, 



where g(x) is an arbitrary function, and q is a two-dimensional vector. Let us substitute 
this relation to (SB) and make the change of variables 



V n + r 2 

Defining an impact parameter p by a relation 



q 



M"ir 2 
n + r 2 



9. 



r x [ri x r 2 ] 



where r = r 2 — ri, taking into account, that at 9 <C 1 the impact parameter p « 0rir 2 /(ri + 
r 2 ) and p r, we can rewrite (p6|) as 



jD (o) ( 



zc 



1*2, ri e 



Aix 2 r 



J dqexp iq 2 — zAr ^ cfeV ^ri + xr + qy / 2r 1 r 2 /fi;r 



(2* 



Here q is a two-dimensional vector perpendicular to vector r. This formula was obtained for 
small angles 9. However, it is correct also for 6 ~ 1, since in this case the term, proportional 
to q the in argument of a potential can be neglected and the Green function (|28|) turns into 
the eikonal one. In particular, for 9' <^ 1 it coincides with (27). The expression ( ^8|) agrees 
with the main (without the correction) term of (|24]) since in the main approximation it is 

possible to replace \Jx{l — x) by y / r 1 r 2 /r 2 . 



3 Wave functions in the quasiclassical approximation 

The obtained expressions for the quasiclassical Green function allow us to find the quasi- 
classical wave functions with the first corrections. The quasiclassical wave functions were 
found earlier in papers || |§ for the case of Coulomb field (Furry-Sommerfeld-Maue wave 
function) and in [M for arbitrary decreasing central potential. These wave functions were 
calculated in the main approximation. In paper jll]| the wave functions and corrections to 



them were found for arbitrary potential in the eikonal approximation. 

To calculate the wave functions, we use the known (see, for example, [I2|) relations 



lim D( Wik) = -^^fe^rO , (29) 
r2^oo 47rr 2 

lim D^\r 2 , ri \e) = - e -^^^\r 2 ). 

ri— >oo 4777"! 

Here pi = —nri/r 1 , p 2 = nr 2 /r 2 , ^ + )(r) ( ^^(r)) denotes a solution of the Klein- 
Gordon equation containing at the infinity a plane wave with momentum p and a diverging 
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(converging) spherical wave. From (|2ll) and (^9|) we obtain for the quasiclassical wave 
function of the Klein-Gordon equation with the first correction taken into account 



(±)/ 



f d(\ f°° 
± / — exp ipr ± iq 2 =F iX / cfeWr,,.) 
J in I Jo 



(30) 



OO X 



x 1 + — y(r) ± l -JdxJdy[^-y] (V ± V(r x ))(V ±V(r y )) 



o o 



r x = r =p px/K + qy 2x/ k , k = |p| . 

Here q is a two-dimensional vector perpendicular to p, Vi is the component of a gradient, 
perpendicular to p. Similarly, for the Dirac equation the quasiclassical wave function with 
the first correction is obtained from (|25[): 



± / — exp ipr ± iq 2 =F i\ [ dxV(r x ) { 1 =F — / rfa; ck W(r x ) + 
J in Y Jo J 2k j 



+ — V(r) ± l -JdxJdy{^-y) (V ± V(r x ))(V ± V(r y )) u x p 



(31) 







where is a free positive-energy (A = 1) and negative-energy (A = —1) Dirac spinors with 
momentum p. In the case, when V(r x ) can be expanded in q, the formula ( |3~TD turns into 
the wave function in the eikonal approximation with correction being in accordance with 



the result of 111 



4 Scattering of charged particles 



In this section we use the obtained quasiclassical Green function for the calculation of the 
small-angle scattering amplitude of the high-energy charged particle in localized potential. 
In the case of a particle with spin zero the scattering amplitude has the form: 



/(p 2 ,Pi) = -£J dr^H*(r)0(r)e^ r . 
Recollect, that = XV - V 2 /2k. 



(32) 



For the sake convenience we put the origin of the reference frame at the point of potential 
maximum and pass to cylindrical coordinates with axis z along the vector p 2 . Substituting 
(BO) in (|32"D, we obtain with the account of the first correction 



IK 

2^2 



dz I dp I (iqexp 



-iQ z z — iQ ± p + iq — iX dxV(x + z,p x 



x (33) 



xXV(z,p) 



CO X 



l+ l -jdxjdy [y/xy - y] (V p V{x + z, p x )){V p V{y + z, p y )) 







where p x = p + J2x/k q, Q = p 2 - p v 



S 



We demonstrate first that in the scattering problem not the eikonal wave function but 
the quasiclassical one should be used. For this purpose we calculate the amplitude in the 
main approximation, that corresponds to the replacement of the factor in square brackets 
of (|33| ) by 1. We split the region of integration over z into two: (—00, 0) and (0, 00). In the 
region from zero to infinity it is possible to neglect quantum fluctuations, which means the 
replacement p x — > p and then to integrate over q. Finally, the contribution of this region 
reads 



/+ = — — / dz dpexp —iQ z z — iQ ± p — iX dxV(x + z, p) 
2tt J J L ~ Jo 



XV(z,p) 



(34) 



Now, integrating by parts over z with the help of the relation 

. d 



-iX 



roo 

/ dxV(x + z, p) 
Jo 



- l d~z 1 6XP 



-iA 



poo 

/ dxVix + z, p) 
Jo 



- 1- 
- 1 



(35) 



(36) 



XV(z, p) exp 
we obtain 

f+ = -7^ J dpexp[-iQ ± p] [exp -iX J dxV(x,p) 

00 

—iQz J dz exp[— iQ z z] ^exp — iX J dxV(x,p) 



The main contribution to the integral over z in this formula comes from z ~ p. The relative 
magnitude of the contribution proportional to this integral, in comparison with the first 
term is Q z p ~ Q z /Q± "C 1, therefore in the main approximation this contribution can be 
neglected. In the region from — 00 up to zero in the main approximation we can replace x 



with \z\ in p x . Performing the shift p — > p — qJ2\z\/K we obtain 



/- 



IK 
2^ 



-iX 



dz J dp J dq exp 

-00 

/ dxVix + z, p 
Jo 



-i(Q z - Ql/2/c) z - iQ ± p + i(q + Q ± J\z\/2kY- 



XV(z,p-qJ2\z\/K) . 



(37) 



Note, that at up ^> 1 at any z the condition J2\z\/k max(|^|,p) holds, which allows us 
to neglect the term proportional to q in the argument of the potential. In the small-angle 
approximation Q z = that is, the term oc z in the exponent vanishes. Taking the 

integral over q, and then over z with the use of the identity (ff5D, we find 

/_ = — — / dpexp[—iQ ± p] fexp — iX J dxV(x,p) —exp — iX J dxV(x,p) j . 

(38) 

Combining / + and /_, we obtain the expression known as the scattering amplitude in the 
eikonal approximation: 



f = — — J dpexp[—iQ ± p] fexp — iX J dxV(x,p) 



(39) 



Usually this result is derived with the use of the eikonal wave function in the whole region 
of integration over z and neglecting the term Q z z in the exponent. It follows from our 
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consideration that at arbitrary momentum transfers both these approximations, generally 
speaking, are incorrect. For example, consider a screened Coulomb potential with the 
radius of screening r c . The main contribution to the amplitude comes from the region 
z ~ r c , p ~ l/Q±. If the quantity Q z z ~ Q\r c /K is not small in this region as compared 

to unity, the term with q in p x can not be neglected, since p ~ 1/ Q± < Jrjn. Thus, the 
eikonal wave function becomes inapplicable in the region of the main contribution to the 
amplitude. Additionally, the term Q z z in the exponent can not be neglected. Keeping this 
term and still using the eikonal wave function leads under condition of Q\r c j /« > 1 to wrong 
result for the scattering amplitude. In particular, acting so, one can not reproduce a well 
known result ( Rat herf or d formula) in the limit r c — > oo (unscreened Coulomb potential). 
Thus, the formula fl39|) is valid for any Q± <C k, however, its correct derivation can not be 
done with the use of the eikonal wave function. 

Let us pass now to the calculation of scattering amplitude with the first correction. For 
the sake convenience we first set the lower limit of integration over z to —L , and then take 
the limit L — > oo. Using the identity 



XV(z, p) + i-^- + A / dx-^= V p V(z + x, p x ) 
oz J J2kx 



exp 


r 

-iX 




J 







(40) 

it is possible to perform the integration by parts over z in the first term of square brackets 
of ([33]). As a result, the expression for / acquires the form f — fo + fi, where 



fa 
fx 



IK 

2^' 



K f f Q I f°° 

— lim dp mqexp i— —L — iQ^p + iq 2 — i\ / dxV(x, p x+L 

7T 2 L^oo J J 2k " J-L 

- Jdz Jdp Jdqexp —i-^z — iQ^p + iq 2 — i\ J dxV(x + z,p x 



(41) 



x 



x 



+ 



^V(z, P) J dx J dy {Jxy - y) (V p V(x + z, p x ))(y p V(y + z, p y )) + (42) 



oo 

A / dx , V p V(z + x,p x 



Qi 

2k 



V2 



The term independent of potential and vanishing at Q ± 7^ 0, is omitted in expression for 
fo- We will take it into account explicitly in the final expression for the amplitude. In 



order to find the limit L — ► 00 in function fo, we make the shifts p — > p — qJ2L/n and 



q — > q — Q±y L/2k. After that the calculation of this limit and integration over q become 
elementary. With the correction of the order of Q±/k taken into account we obtain 



fo 



iK r f°° ia r 

— — / dpexp —iQ ± p — i\ / dxV(x,p) 1 + — / dx xQ ± 'V p V(x, p 
2n J ~ j-00 2k J 



(43) 



The integral over p in the expression for f\ (|42|) converges at p ~ l/Q±. The quantum 
fluctuations are important only at z < and x close to —z. Since the quantity j\ is small 
being a part of the correction, p x in the formula (f42|) can be replaced by p^i. Besides 



the factor V(z, p) can be replaced by V(z, p\ z \) because the difference between these two is 
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small at any z. The integral over z from zero to infinity converges at z < p, so zQ, 2 l /2k ~ 
pQ, 2 ± /n ~ Q_l/k <C 1. Therefore, it is possible to replace in the exponent zQ 2 ± /2k by 
— \z\Q, 2 ± /2k. After integrating by parts over q the term proportional to q this vector enters 
everywhere, except the exponent, only in the argument of potential. Then the expression 
for fi can be transformed to the form 



fx = J dz J dp J ^ exp 



' 2k 



\z\ - iQ ± p + iq 2 



9(z,p+ J2\z\/nq 



(44) 



with some function g. If we perform the shift p — > p— qy 2\z\/k the integral over q becomes 
elementary. As a result we have 



oo 

h = J dz J dpexp[-iQ ± p]g(z,p) ■ 



(45) 



Finally, we obtain the following expression for f\ 

oo 

1 



fi = / dp [ dzexp —iQ, ± p — i\ [ dxV(x + z, p) 

2n J J L ~~ Jo 



x 



(46) 



x 



+ 



i\V(z, P)jdxj dy(y/xy-y) (V p V{x + z, p))(V p V(y + z } p)) 
o o 

oo oo 



When deriving this formula we have integrated by parts over p the term proportional to 
Ql- As can be checked, the integrand in (|46|) is a total derivative over z and the integration 
over z becomes trivial. Combining obtained expression with (|43|) , and integrating by parts 
over p the term proportional to Qj_, we obtain the scattering amplitude with the first 
correction 



%K f \ /"oo 

/ = — — / dp exp [— iQj_p] { exp — iX / dxV(x,p) 

2n J [ L J-oo 

n r , °° 

1 A c 

dxV(x,p) — / dxxA p V(x 1 p) — 

3 J 2k J 



(47) 



+ exp 



K 



X 

dx J dyy(V p V(x,p))(V p V(y,p)) 



Here A p = Using the wave function (|3l|) , it is possible to show (see, for example, 
|T3||), that the small-angle scattering amplitude for particles with spin 1/2 coincides with 
amplitude ( f47D for spin- zero particles including the terms of the order of Qj_/k. 



The obtained amplitude has correct properties with respect to shifts. It follows from 
(j33D, that after the replacement V(r) — > V(r + r ) the amplitude /(p 2 ,p 1 ) acquires the 
factor exp[zQr ]. The amplitude (|47| ) obviously has this property with respect to the shift 
in the direction, perpendicular to p 2 , i.e., to z-axis. Let us consider now the shift along 
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the z-axis: V(z,p) — > V(z + z ,p). Performing the change of variables X — > X — Zq and 
y — > y — z hi the integrals over x and y , we obtain the additional term proportional to z : 



/ (ipexp —iQ ± p — i\ / dxV(x,p 



CO 

oo 



X 



(4* 



y A p V A (x, p) — % | y <ixV p V(x, p) 



47T 



dpexp 



— «Qj_p] A p I exp[— zA / dxV(x,p) 

\ J — oo 



- 1 



Integrating by parts over p and substituting Q 2 l /2k — > Q z , we hnd for Q^^o *C 1 that 
/ + 5/ ~ / exp(iQ 2 -2o)- Thus, the correct transformation properties of the scattering 
amplitude with respect to shifts holds with the same accuracy, as formula ([47]) itself. 



For the potential satisfying the condition V(z, p) = V(—z, p) the expression (|47D agrees 
with that obtained in [O] in the eikonal approximation with the account for the correction. 
However, as explained above, the formula ( 47|) holds even when the eikonal approximation 
is inapplicable. 



5 Delbruck scattering 

The process of coherent photon scattering in the electric field of atoms via virtual electron 
-positron pairs (Delbruck scattering) has been intensively investigated both theoretically, 
and experimentally (see, e.g., the review |14j])- I n this section we consider the forward 
Delbruck scattering in a screened Coulomb field as one more example of the application of 
the quasiclassical Green function obtained. The photon energy uj is assumed to be large in 
comparison with the electron mass m. According to the optical theorem, the imaginary part 
of this amplitude is proportional to the cross section of electron-positron pair production 
by a photon in the field of atom. To describe the propagation of light in matter it is also 
necessary to know the real part of this amplitude. 

For the calculation of amplitudes of Delbruck scattering with the help of Green functions 
one is forced to use different approximations for these functions for different momentum 
transfers A = k 2 — \<L\ (k 1; k 2 are the momenta of initial and final photons). At A ~ u the 
quasiclassical approximation is inapplicable, since the basic contribution to the amplitude 
is given by the orbital moment I ~ uj/A ~ 1. For such momentum transfers the amplitude 
was calculated in |15 |. In region uj ^> A ^> m 2 /uj it is possible to use the quasiclassical 
Green function in the main approximation |16[ . The calculation of Delbruck scattering 
amplitude in Coulomb field at A < m 2 /u required a special consideration (see |17|, H). 



In this case the contribution to the amplitude is given by the impact parameters p up to 
p ~ uj/m 2 . At such impact parameters it is necessary to take into account a correction to 
the quasiclassical Green function 0. For the screened Coulomb potential with uo/m 2 ^> r c 
( r c is the radius of screening, in Thomas-Fermi model r c ~ (Ma)~ 1 Z~ 1 / 3 ^> 1/m) the 
contribution to the amplitude is given by the impact parameters p < r c uj/m 2 . In 
this case, the quasiclassical Green function without corrections can be used for arbitrary 
A <C co [fl [/J. However, if A, r^ 1 < m 2 /u), the correction is very important. Moreover, 
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the expression for the forward scattering amplitude (A = 0), obtained with the use of the 
quasiclassical Green function without corrections is, strictly speaking, indefinite. Below 
we derive the amplitude of forward Delbruck scattering at arbitrary ratio between r c and 
uj/m 2 . 

As shown in ||, the amplitude of forward Delbruck scattering for high photon energy 
can be presented as follows 



M = %ol J de J dvidv2 exp[zk(ri — r 2 )] x 



(49) 



xSp 



[2e*p 2 - e*k)D(r 2 , Ti\u - e) 



(2epi + ek)D(Ti, r 2 | — e) 



where e, k are polarization and 4-momentum vectors of a photon, pi 2 = — zVi )2 - In this 
formula the subtraction from integrand of its value at zero external field is assumed to be 
done. Since in the case of central field the amplitude of forward scattering does not depend 
on polarization of photon, it is convenient to make the substitution e*ej — > (5ij — kikj/u 2 )/2 
in ©. 



Let us pass in (HS) from variables r 12 to 



r2 - r 1; 



r x [ri x r 2 ] 



Since pr = 0, the integration over p is carried out in a plane, perpendicular to r. The 
main contribution to the amplitude comes from the integration region, where r ~ u/m 2 , 
\z\ ~ 1 and angles between vector r and k are of the order 5 r ~ m/w « 1. Due to the 
smallness of angles 9 r it is possible to consider vector p to be perpendicular to k too. 
Besides, it is obvious that the main contribution is given by p < min(l, m 2 r c /uj). 

Let us split the region of integration over p into two: from up to po and from po to 
oo, where m/u <C po min(l, m 2 r c /uj). In the first region (at p < po) the following form 
of the quasiclassical Green function can be used 



-D(r2,ri|e) 



le 



4ir 2 r 
x exp 



iq 



1 + 



a. q 



KT 



i\r 



e y 2r\r 2 
J dxV (y\ + xr + qy / 2rir 2 /fi;r 



(50) 



which we obtained by substituting (28) in (y) and performing some transformations of 
the term containing a -matrix. Namely, the longitudinal components of gradient, which 
in comparison with transverse ones have additional smallness p, were omitted and the 
integration by parts over q was performed. Besides, by virtue of the definition of po, in this 
region we can neglect screening and replace V(r) by a Coulomb potential V c (r) = —Za/r. 
After this the integral over x in the exponent can be easily taken. 

The screening is essential only in the second region, where we can use the representa- 
tion (|13|) , that is, the eikonal Green function with the first correction. The substitution of 
(|13|) into fl49D results in the cancellation of the potential-dependent terms in the exponent. 
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Therefore, the contribution from the second region does not contain any powers of the 
potential except the second one, which corresponds to the first Born approximation. Be- 
sides, due to the phase cancellation we must take into account the correction to the Green 
function. In this region the contributions to the amplitude rising from the correction and 
from the main term in Green function turn out to be of the same order of magnitude. 

Going over to the calculation of the contribution Mi from the first region, we substitute 
the Green function (|50|) for a pure Coulomb potential into (49) , differentiate and take 
the trace. Using the smallness of angles between the vectors r and k and expanding 
\/e 2 — m 2 ~ \e\ — m 2 /2\e\, we obtain 



Mi 



deeK drr dO r 




dp 



dz 



P<P0 



{z{l-z)f 



Jd(\id(\2 x 



(51) 



2iZc 



exp < i- 



uOl- 



m 2 uj 



en 



+ 



z(l-z) 



X 



X 



2e«;[qiq2 - z(l - z)Ol) + 



CO 



■(eqi - Kq 2 , qi 



q2j - 1- 

r 



4(s(l -z)) 

where k = uo — e as well as vectors qi i2 , and r are two-dimensional vectors, perpendicular 
to k. Note that the integral over z in this formula is taken in the limits from zero to unity. 
The reason is that outside this interval the Green function have the eikonal form and the 
phases depending on a potential are cancelled. For p < p <C 1 it results in a negligible 
contribution from the region outside the interval < z < 1 compared to the contribution 
from this interval. 



Let us integrate now over 6 r , pass from the variables q i)2 to Q = (qi + q 2 )/2 and 
q = (qi — q 2 )/2 and make the shift p — ► p + Q. After that the integral over p acquires the 
form 



J 



dp 



Re 



2iZa 



p + q 



i 



(52) 



ip+Qi<po 

The main contribution to the amplitude ( ]5"TD comes from the region Q, q ~ m/uj p > 
where we can neglect Q in the limit of integration in (|52|). To take this integral (^) we 
subtract and add to the integrand the function — 2(Za) 2 [2pq/(p 2 + q 2 )] 2 , which is easily 
integrated: 



Ji 



dp 



P<PO 



-2(Za) 



2pq 

n 2 + q 2 



-4n(Za)Y ln^-1 
V Q 



(53) 



In turn, dealing with the difference we can extend the integral over p to infinity owing to 
the fast convergence of the integral. To calculate this integral, it is convenient to multiply 
the integrand by 



1 a A 2 /°q 

-l V p z + q z 



J J -i \y\ 



(W - 1)) , (54) 



and change the order of integration over p and y. Integrating over p, we obtain 



J 2 = 4nq 2 



dy 

y 3 



Re 



w 

1 + 3/, 



iZa 



- 1 + 2(Za) 2 y 2 



(55) 
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Using the replacement y = tanh r we find for J = J\ + J 2 



J = 8nq 2 (Za) 



Po 



1 + Re<0(l + iZa) + C 



(56) 



where C = 0.577... is the Euler constant, ip{x) = d\nT(x)/dx. It is convenient to take the 
remaining integrals in the following order: over Q, q, r, z, and e. Finally, the contribution 
from the first region reads 



Mi 



.28a(Za) 2 uo 
9m 2 



In 



ujpo 



m 



7T 



i— — Re0(l + iZa) 



4T 

C 

42 



(57) 



The contribution of the higher orders of the perturbation theory in the external field 
(Coulomb corrections) is given by the term — Re^(l + iZa) — C in (|57|) and coincides 
with the known result [|Tj]]. Thus, the Coulomb corrections are completely determined by 
the first region, when the quasiclassical Green function is not reduced to the eikonal one. 

Let us pass to the calculation of the contribution M 2 from the second region. Taking 
the derivatives over r 12 , calculating the trace over gamma-matrices and integrating over 
6 r we come to the following representation for 



Mo 



a 



2tvuj 



3 

de Jdrr 2 exp 



ujrm 

i 

2sk 



1 1 



dp I dz dx dy x 



(5? 



2y(l - x)[2#(x - y) + 1] 



J p>po 

" 2en 







[V P V(R Z - X )} ■ [V P V(R Z . 



yji ) 



where R s = r\/s' 2 + p' 2 . In this formula the terms, antisymmetric with respect to the substi- 
tution e ^ uj — e, z ^ 1 — z , are omitted, since their contribution to the integral vanishes. 
We emphasize, that in contrast to the first region, in the second one the integration over 
z is performed in infinite limits. Changing the variables z^z + x, y^-y + x in triple 
integral over x, y, z and performing the integration over x, we have 



M, 



a 
2t\uj 



de drr 2 exp 



x 







ujrm 



—i- 



UJ 



2en 



J dp J dz J dy(l -y)x 



(59) 



P>P0 



y) 2 + 2y- — 

SK 



[V P V{R Z )\ ■ [V p V(R z . y )} . 



At r r c the potential V(r) w —Za/r. Therefore, it is convenient to present the contribu- 
tion (H) as a sum M 2 = M 2 C + 5M, where M 2 C is the value of M 2 at V(r) = V c {r) = -Za/r. 



At V = V c the integrals over r and e can be easily taken, and we find 



M 9 C 



2ia(Za) 2 uj 



dpp J dz J dy(l - y)- 

PO —oo 

Using the Feynman parametrization of denominators 

„ i 



2(l-y) 2 /9 + y/3-l 
z 2 + p 2 f/ 2 [{z- y) 2 + p 2 } 3 / 2 ' 



(60) 



(AB) 3 /2 



7T 



v(l 







[Av + B(l - v)Y 
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we take the integrals over p and z, and then over y and v bearing in mind, that po ^ 1- 
Finally, for this (Coulomb) contribution we obtain 



M 2 C 



.28a(Za) 2 u ( p 31 



9m 2 



In — H 

2 21 



(61) 



The sum of the contributions fl5"?D and (|6lD gives the known result for a pure Coulomb 
potential jp!7[| : 



,28a Zafw [, 2w ,tt „ , , _ 109 

M c = i ^— - — In i Reib(l + iZa) -C 

9m 2 I m 2 ^ v ; 42 



(62) 



In the term 8M, connected to screening, the lower limit of integration over p can be 
replaced by zero owing to the convergence of the integral at p — > 0. Using the momentum 
representation for potentials 



V{r) 



dp 



■e tpr V(p) 



we take the integrals over p and z. The result of this integration is proportional to 8(jp\ 
P2). Integrating over P2 and over the angles of vector pi, we find 



5M 



2tt 3 
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de drrexp 



—1- 



2ek 



x / dp v V\v) - (4vrZft) z 



dy{l - y) 
( sin ( cos C 



|(l-Z/) 2 + 2y-- 



x (63) 



V C 3 



c 2 



where ( = rpy. Passing from the variable r to (, integrating by parts over p and taking 
the integral over (, we obtain 



8M 



auo 
27T 3 m 2 



1 1 
dxx(l — x) \ dy(l/y—l) 



4, „ 1 

(l-yf + 2y-— 

3 a;(l — x) 



(64) 



x / dp(d pP 4 V 2 (p)) 



- + - 1 - - In 



1 +7/ 



?7 



1 — 77 — z0 x 



where the substitution e — > cux is made, 77 = 2ux(l — x)py/m 2 . Now we make the change 
of variables y — > y/(2x(l — x)), change the order of integration over x and y, and take the 
integral over x. Finally, the contribution to the amplitude of forward Delbruck scattering 
due to screening reads 



5M 



1/2 



187r 3 m 2 



dp{d pP *V\p)) / *H 



- + - 1 - - In 



1+7/ 



1 — 77 — i0 



(62/ + 7y + 7) Jl - 2y + 3y(2y 2 -3y-3) In 



1 + 

1- y/T^, 



(65) 
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where rj = uopy/m 2 . This formula holds for arbitrary form of a screened Coulomb po- 
tential. In some special cases it can be essentially simplified. For the potential V(r) = 
— Zaexp(— f3r)/r, when V{p) = —AmZa/{jp 2 + (3 2 ), all integrals in (|65|) can be taken: 



SM 



Aia(Za) 2 u 
9m 2 



33 - 13r 2 + -r 4 + -r(24 - 13r 2 + 3r 4 )L 



'1 + 



2im 2 



In 



T- 1 
T+ 1 



For more realistic Moliere potential [ITH], we have 



9t 2 + t 6 )L 2 



(66) 



V(p) = -AnZa ^ 



a\ = 0.1, a 2 = 0.55, a 3 = 0.35, (3 n = f3 b n , 
b x = 6, b 2 = 1.2, 6 3 = 0.3, fa = mZ 1/3 /121. 



(67) 



In this case 
6M 



1/2 

4ia(Za) 2 u f dy 

y 



9m' L 



(6y 2 + 7y + 7)Jl-2y + 3y(2y 2 - 3y - 3) x 



(68) 



x In 



l + y/1^2^ 
1- V^%J 

where 7„ = u)f3 n y/m 2 . 



1 - 



2?; 



2 . (Jn + i 
In 



7n + 7fc 7n-7fe \7fc + ^ 



In + 1 



As known, the imaginary part of the amplitude of forward Delbruck scattering is con- 
nected with the total cross section a of the electron-positron pair production by photon in 
an external field via the relation a = ImM/u. For the Moliere potential, using (|62|) and 
(|68|). we obtain the following expression for o: 



1/2 



dy 

y 



a 



,6 



28a(Za) 2 f 2co 



9m 2 



109 

In — - Reip(l + iZa) -C 

m 42 



^y 2 + y + 1)^1 -2y + -y(2y 2 - 3y - 3) In 



1- y/T^Zy-, 



x 



X 



ll-lk 



In 



1 + 7^ 
,l + 7 fc 2 , 



7; 



+ 7i; 



(69) 



which agrees with the results, known in literature (see, e.g. [|T^ 



Let us discuss now the dependence of the real part of Delbruck amplitude from the 
photon energy 10. When uo/m 2 <C r c the screening can be neglected and ReM = ReM c = 
147ra(Za) 2 LV /9m 2 . The linear growth of the real part with increasing uo gradually becomes 
weaker and for uo/m 2 3> r c we obtain from ( |6^) and (|65D 



ReM 



a 2 



2vr 3 



CO 



m 2 r, 



dp 
P 



(d pP 4 V 2 (p)). 



(70) 
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Here only the term containing a higher degree of the large logarithm is kept. For the 
illustration the real part of the forward Delbruck amplitude is shown in Fig. 1 for Z = 82 
and V(r) = —Zaexp(—maZ 1 / 3 r)/r as a function of oj. 

One of the basic mechanisms of elastic scattering of a photon is the Compton scattering 
on atomic electrons. For the forward scattering the amplitude of this process is real and 
does not depend on uj: Mcomp = —^nZa/m. It is seen in Fig. 1, that the interference 
between the amplitudes of Compton and Delbruck forward scattering should be taken into 
account already at relatively small energies. 

From the calculation of the Delbruck scattering amplitude we learn ones more, that the 
use of the eikonal approximation for the description of high-energy small-angle scattering 
processes without proper ground can lead to incorrect results. For instance, the Coulomb 
corrections in the imaginary part of the forward Delbruck scattering amplitude (and, there- 
fore, in the total cross section of pair production) would be completely lost if we used the 
eikonal Green function in the calculation. 
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1. The real part of the forward Delbruck scattering amplitude for the potential 
— Za exp(— maZ l ^r) /r in units of AitZa/m, Z = 82 as a function of to. 
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Fig. 1. 



